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Abstract 

Pointing out the incompleteness of conventional macroscopic Maxwell equations (M-eqs.), 
we propose a new form derived from the long wavelength approximation (LWA) of micro- 
scopic nonlocal response. From the general Hamilonian of matter and matter-EM field 
interaction (containing spin dependent terms due to relativistic correction), we first set up 
the simultaneous equations for microscopic "vector potential A and induced current den- 
sity /" , and then extract the macroscopic components by applying LWA. This leads to new 
macroscopic M-eqs. with a single macroscopic susceptibility Xcmik, lu) between I and A, 
which describes both electric and magnetic polarizations and their mutual interference in its 
fully quantum mechanical expression. In the absence of chirality and under the condition to 
use magnetic susceptibility defined with respect to B, this scheme is shown to be equivalent 
to the conventional " e and fi scheme" . In the case of chiral symmetry, the phenomenological 
constitutive equations by Drude, Born and Fedorov cannot be justified by this microscopic 
approach. As a single susceptibility scheme of macroscopic M-eqs., this result is on an ad- 
vanced level by its fully quantum mechanical description of the whole set of 0(fc°), O(fc^) 
and O(fc^) terms of susceptibility providing a consistent picture of such a scheme. 

1 Introduction 

Maxwell equations (M-eqs.) are one of the most fundamental theoretical frameworks of physics, 
dealing with the interaction between electromagnetic (EM) field and charged particles with 
magnetic moment, which describes the EM response in both static and dynamical regimes. 
They consist of four laws, i.e., Gauss laws for electricity and magnetism. Ampere's law and 
Faraday's law. There are two forms of them, i.e., microscopic and macroscopic M-eqs., and both 
are actively used for research and education. The fundamental formulas and definitions about 
micro- and macroscopic M-eqs. are given with usual notations in Appendix A. 

Historically, macroscopic M-eqs. appeared first as a consistent set of the four equations 
mentioned-above for macroscopic matter systems. It was a phenomenology, since nothing was 
known about electrons, quantum mechanics and relativity at that time. As the particle picture 
of matter was established, one could correctly transform macroscopic M-eqs. to microscopic 
ones describing the EM field in microscopic systems [Ij. The microscopic form, combined with 
quantum mechanics and relativity, lead to quantum electrodynamics (QED), which is the most 
accurate theory of physics in precise agreement with experiments [21 13] . 

On the other hand, the particle picture of matter in combination with quantum mechanics has 
been used to promote macroscopic M-eqs. from a phenomenological level to a semiquantitative 
one by treating susceptibilities, not as parameters, but as the quantities to be calculated quantum 
mechanically. It is based on the viewpoint that macroscopic behavior of EM field must somehow 
be derived from the averaged motions of microscopic EM fields and charged particles. This type 
of rewriting has been called the derivation of macroscopic from microscopic M-eqs., and such 
semiquantitative M-eqs. in terms of the field variables "E, D, B, H" have been used with much 
success for the studies of macroscopic EM response in a vast area. 
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However, in spite of the long history with great success, the macroscopic scheme in terms of 
"E,D,B,H" seems to be incomplete in comparison with microscopic one with respect to the 
follwing points. 

1) Lack of explanation why the number of susceptibility tensors are different betweeen mi- 
croscopic and macroscopic M-eqs. 

2) Lack of uniqueness in dividing charge density into true and polarization charge densities 

3) Lack of uniqueness in dividing transverse current density into the contributions of magnetic 
and electric polarizations 

4) Different wave number dependence of magnetic permeability fj, for two types of magnetic 
dipole (Ml) transitions, i.e., spin and orbital Ml resonances 

The details of each point are explained below. 

1) In the derivation of macroscopic M-eqs. from microscopic ones, which can be found in 
many textbooks of electromagnetism [4], it is admitted from the beginning that the two 
vector fields, i.e., electric and magnetic polarizatoins, are the fundamental quantities of 
macroscopically averaged matter. This is, however, not obvious, because in the more basic 
microscopic M-eqs. the dynamical variable of matter for EM response is a single vector 
field, i.e., "current density". (Charge density is not an independent variable, since it is 
related with current density via continuity equation.) This means that the number of 
constitutive equations (or susceptibilities) is different between the microscopic and macro- 
scopic M-eqs.. It is not obvious that macroscopic averaging requires such an increase in 
the number of constitutive equations [5j . 

For systems with chiral symmetry, where elecrtic field can induce elecrtic, as well as mag- 
netic polarization, (and magnetic field can induce magnetic, as well as electric polariza- 
tions,) a phenomenological approach by Drude-Born-Fedorov (DBF) constitutive equations 
is known to be useful [6], [7]. (See Appendix A.) This scheme requires an additional (chi- 
ral) susceptibility. Then, the number of suseptibility in macroscopic M-eqs. is increased 
to three, in further contrast with microscopic M-eqs. 

2) It is obvious that the separation of a given charge density p{r, t) into true and polarization 
charge densities, pt{r, t) and Pp{r, t), respectively, cannot be done uniquely. Since the only 
constraint is the charge neutrality J dr pp = 0, there are infinitely many ways to rewrite 
p into pt + Pp. Because of the continuity equation, this problem can be rephrased as the 
non-uniqueness of separating the longitudinal (L) component of current density J into 
the conduction current density Jc and the L component of the current density dPi^/dt 
induced by electric polarization P. 

3) A similar problem exists in separating the transverse (T) component of J into magnetiza- 
tion (iW) induced current density c V x M and the T component of i-'-induced current 
density dP'^/dt. In order to do this division uniquely, we need to define P and M pre- 
cisely, which, however, is not always possible. Beause in chiral symmetry the electric 
dipole (El) and magnetic dipole (Ml) transitions are mixed, which allows P induced by 
magnetic, as well as electric field and M induced by electric, as well as magnetic field. 
Then, there is no unique way to separate P and M. 
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4) The problem about ;^(u;) is the absence of its general expression including the contributions 
from both spin resonance and orbital magnetic dipole (Ml) transitions. Since the deviation 
of n{io) from unity is due to the Ml transitions of matter system, one would expect similar 
forms of resonant terms for both transitions. However, the conventional ways to write the 
two contributions are different with respect to their dependence on the wave number of 
the resonant EM waves. In the case of spin resonance [8j, a typical magnetic susceptibility 
has the form Xm = b/[uJo — uj — ij], where cuq is the transition energy of a spin resonance, 
and b the squared transition matrix element of spin magnetization operator. The resultant 
expression /x = l+47rxm does not contain the wave number of the resonant EM field except 
for uj in the denominator. On the other hand. Ml component of orbital transitions, which 
occur in any energy regions including nuclear transitions (e.g., Mossbauer transitions), 
produces an EM field with an amplitude proportional to the product of the resonant 
wave number and the magnetic dipole moments of the transition, according to the 
standard multipole expansion of the EM field [9]. This means that the /u(u;) in this case 
has a resonant term of intensity proportional to k'^. Thus, the two cases of Ml transition 
apparently show different k^ dependences of /x(a;). 

Another difficulty about /x is that, in the standard dispersion equation {ck)'^ /uP = e/i, 
there appears the product, not the sum, of electric and magnetic contributions in spite 
of linear response treatment. If we consider the possibility of the mixing of electric and 
magnetic excitations in chiral symmetry, there is no essential distinction between them as 
matter excitations. Therefore, it would be more reasonable to have their contributions as 
a sum, rather than as a product, in the dispersion relation of linear response. Though 
this problem looks different from the one in the previous paragraph, it turns out later that 
they have a common origin. 

To the author's understanding, the origin of the incompleteness mentioned above lies in the 
use of P and M as the independent, fundamental variables in deriving macroscopic M-eqs. This 
independence is questionable, if we think of chiral symmetry where electric and magntic transi- 
tions contributing to P and M are mixed with each other. As a logic to derive macrosacopic 
M-eqs. in the conventional manner, it is necessary to show the general definition of P and M 
from the first principles, i.e., from the precisely defined matter Hamiltonian and matter-EM 
field interaction. As will be discussed in Sec. 3.1, one would realize the incompleteness from the 
very effort to derive the matter-EM field interaction in an appropriate form to calculate electric 
and magnetic susceptibilities from the ffi'st principles. 

Therefore, in order to find right answers to the problems mentioned above, we need to use a 
new and more reliable logic in deriving macroscopic M-eqs. As such, we start from a microscopic 
response theory, a most reliable basis for the EM response of matter in semiclassical regime. 
In this scheme, all the details of matter states are described quantum mechanically, and all 
the wavelength components of EM field are taken into account in a semiclassical way. Matter 
Hamiltonian and matter-EM field interaction are precisely defined from the Lagrangian for a 
general interacting system of charged particles and EM field, and explicit spin dependence is 
considered as relativistic corrections. In the microscopic response theory, we solve, so to speak, 
Schrodinger eq. and micoscopic M-eqs. selfconsistently, or more specifically, the coupled integral 
equations for induced current density and vector potential, which contain all the wavelength 
components. This way of solution becomes feasible due to the separability of the microscopic 
nonlocal susceptibility as an integral kernel, as will be shown in Sec. 2. We call this scheme 
"microscopic nonlocal response". 

Then, as a straighforward procedure of macroscopic averaging, we apply long wavelength 
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approximation (LWA) to the coupled integral equations. Thereby, we do not assume or request 
to recover the conventional form of macroscopic M-eqs. but just take the outcome as a logi- 
cal consequence of this approach and compare the result with conventional form. Here also, 
the separability of microscopic susceptibility plays an essential role to give an explicit form of 
macroscopic susceptibility quantum mechanically. This type of derivation has never been done 
before to the author's knowledge, and in fact it gives reasonable answers to all the problems 
mentioned above. It should be noted that this approach describes the macroscopic averaging, 
not just as a concept which might allow various interpretations, but as an unambiguous mathe- 
matical procedure of LWA applied to the well-defined integral equations of microscopic nonlocal 
response. This enables us, in principle, to check the validity condition for the approximation, 
since we clearly know what is discarded by LWA. This point will be discussed in Sec. 3. 5. 

The most remarkable result is that the M-eqs. remain in the same form after applying 
LWA, and the only difference between the micro- and macroscopic M-eqs. lies in the form of 
the susceptibility. Among several proposals of the single susceptibility scheme of macroscopic 
M-eqs. |10 1 lll t ri2]. the present one is on an advanced level in the sense that all the components 
of susceptibility are given in a qunatum mechanical way, which, in principle, allows microscopic 
evaluation, and that it is a logical consequence from the more basic microscopic response scheme. 
Detailed discussions will be given in Sec. 3. 2. 

The existing schemes of M-eqs. in the semiclassical regime are classified into the following 
hierarchy in the order of decreasing accuracy: (i) microscopic nonlocal response, (ii) the new 
macroscopic scheme, (iii) the conventional "e and ^u" scheme. Via LWA, (i) reduces to (ii), and 
(iii) is obtained from (ii) in the absence of chirality and by the use of magnetic susceptibility xb 
defined with respect to B. The justification of the step from (ii) to (iii) goes in a clear cut way as 
shown later, but a similar argument is not possible for the justification of the DBF equations for 
chiral symmetry, so that the DBF eqs. should stay outside the hierarchy as a phenomenology. 

The use of xb is, not only a logical consequence, but also the solution of the problem 4). 
Rewriting the dispersion equation (cfc)^/u;^ = e/i in terms of = 1/(1 — Aitxb), obtained from 
M = XmH = xbB and fi = 1 + 47rxrrn we immediately see the solutions for the both problems 
of 4), which will be discussed in Sec. 2. 4 in more details. The requirement to use xb will be a 
considerable impact on the study of LHM. In Sec. 3. 4 an observable difference is demonstrated 
in the forms of dispersion curve and reflectivity spectrum in a LHM situation for the cases of 
assigning magnetic excitation energy to the pole of Xm and xb- 

Macroscopic M-eqs. are still an important tool in many research fields, such as photonic 
crystals [13], left-handed materials (LHM) or metamaterials |14] . near field optics [15], and so 
on. For this reason and also for its importance in the elementary course of physics education, 
it is worth looking for a more complete form of macroscopic M-eqs. Especially the present 
investigation, which raises a question [5] about the validity of the Veselago's condition for LHM 
"e < 0,// < [16], has direct relevance to the today's very popular subject, LHM. The results 
of this work would affect (at least) the research fields mentioned above on the very fundamental 
level. 

This paper is arranged as follows. In Sec. 2, we give a precise definitin of matter Hamiltonian 
and matter-EM field interaction, set up the simultaneous integral equations of microscopic non- 
local response, and, by applying LWA, extract the macroscopic components, which leads to a 
new macroscopic susceptibility Xem- The dispersion equation in terms of Xem is derived, and the 
conditions to reproduce the conventional one are studied. Discussions are given in Sec. 3 about 
(i) the interaction Hamiltonian in the form explicitly containing E and B, and the condition 
for defining Xc and Xm separately, (ii) the comparison with other macroscopic theories, (iii) the 
comparison with the result of DBF eqs. for chiral symmetry, (iv) the consequence in the LHM 
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studies, (v) the validity condition of LWA, and (vi) simplification of Xcm with parameters. In 
Sec. 4, we summarize the important results. 

2 Formulation 

2.1 Precise definition of matter, EM field and interaction 

We start from the well-known Lagrangian for charged particles in EM field 

1 .,ldA 



C = J2 - ^i<t>{rd + ^v, . Air,) } + J dr ^{(^^ + V^)^ - (V x Af} , (1) 



where A and <^ are the vector and scalar potentials, and ee, ri and V£ are the charge, coordinate 
and velocity of the i-th particle, respectively. From the least action principle for this Lagrangian, 
we obtain the fundamental equations of the system, i.e. the microscopic M-eqs. and the Newton 
equation of each particle under Lorentz force. 

The microscopic M-eqs. are given in Coulomb gauge as (V • A = 0), as 

V^<A = -47rp, 1^-V2A = ^Jt, (2) 

Ot'^ c 

where Jt is the T component of J, i.e., V • Jt = 0, and the charge density p and current density 
J are defined as 

/o(^) = J2^<^ ^ir-ri) , (3) 

e 

•^(f) = ^ ee ve 6{r - ri) . (4) 

e 

The matter Hamiltonian including the interaction with EM field takes the form 

e^rui c ^ i^fr^ \re-re>\ 

in Coulomb gauge, where momentum p is defined by eq.(|57p. The A-linear part of this Hamil- 
tonian 

i^intj = -^ J dr J{r)-Air,t). (6) 

represents the linear light-matter interaction, and the A-independent part is the pure matter 
Hamiltonian. 

If one needs relativistic corrections, we add spin-orbit interaction, mass velocity term, and 
Darwin term to the matter Hamiltonian. Spin Zeeman term can be added to //jntj because of 
their similar form as 



Hsz 



- J dr Mspi„(r) • B{r) = -\ j J.^^r) ■ A{r) , (7) 

where the current density due to spin magnetization is 

Jspm(r-) = cV X Mspin(r) (8) 

and 

Mspin(r) = E Pi^i ^{f - re) , (9) 
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Pe^e being the spin magnetic moment of the i-th particle with spin si. 
Thus, the matter Hamiltonian is 



including the relativistic correction Hj-c, and the matter-EM field interaction H[^t is 

Hint = -l J dr I{r)-A{r,t) , (11) 

where 

I{r) = J{r) + Jspin(r) , (12) 

is the generalized current density operator, containing the contributions from both electric po- 
larization and "spin and orbital" magnetization. The microscopic M-eqs. including this gener- 
alization is 

1 d'^A 2 . 47r 



^'A = -Ir^ (13) 

How to handle static EM field will depend on problems. If we calculate static susceptibility, 
we will put the corresponding term in H\^i. If we are interested, e.g., in electron spin resonance 
or cyclotron resonance, however, we should put the static field part in Hq and the time dependent 
part to -ffint- 

Another point of importance in considering EM response is how to treat L electric field 
induced in matter system. It can arise from true charge density and electric polarization. Let 
us assume for simplicity that matter system has no true charges {J2e^^ — 0)- Its dynamically 
perturbed state is described by oscillating charge and current densities. The polarization charge 
density pp induces L electric field as 

El''''\r) = -vJ dr-'^J^ , (14) 

which is a part of Maxwell field E. Its interaction energy with the polarization of the matter 
can be rewritten as the Coulomb interaction energy among the induced charge density, 

HintL = -J dr P{r) . Et''\r) , (15) 

d./dr-'^li^!:^. (16) 



This leads us to a dichotomy, whether (A) we regard E^^'^'' as an internal field of matter, 
taking full Coulomb interaction into the matter Hamiltonian, or (B) we regard E^^^^ as a 
part of the external source field to induce polarization, taking the Coulomb interaction (jl6|) 
as the interaction between the "external field" and matter polarization, i.e., omitting it from 
the Coulomb interaction energy of the matter. Since the omitted part of the Coulomb energy 
contributes to the energy difference between the T and L mode configurations, the scheme (A) 
contains the different T and L mode energies as the pole positions of susceptibility, while (B) 
reveals the effect of LT splitting, not in the susceptibility, but in response spectra. Though 
the two schemes lead to the same observable result, there is a considerable difference in the 
intermediate steps. 
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In this paper, we take the viewpoint (A). Namely, we define the matter Hamiltonian as the 
sum of kinetic energy and (complete) Coulomb interaction of all the charged particles, and E^^'^^ 
is not considered as external field. Thus the Poisson equation contributes to the internal motion 
of matter, relevant to its quantum mechanical eigenstates, and only the second equation of ^ 
represents the interaction between matter and EM field. (If we were to take the viewpoint (B), 
we would treat -ffjnt + -ffintL as the interaction term, and omit (jl6p from the Coulomb interaction 
of eq.([5]). Using these matter Hamiltonian and matter-EM field interaction, we could calculate 
the time development of the operators /, P, M according to the method in the next section.) 

Also we assume that the incident EM field is T mode, i.e., we do not consider the excitation 
by external charged particle. Thus, in the Coulomb gauge, which is a natural choice for matter 
systems in non-relativistic regime, the matter Hamiltonian contains the full Coulomb interaction, 
and the external EM field is T-mode alone, i.e., A or {-Et andS}. After determining the 
response, we can calculate from the induced charge or current density. (If we were to 
consider the excitation by external charged particle, which will be treated in a later publication, 
we would add an extra term 

= - J dr P{r) • Et'\r) , (17) 

to -ffint- The external L field E^^^^ and an external charge density p(r')ext are similarly related 
as in eq. (ll4p . In terms of the matter-EM field interaction, H[^t + -f^h^L' could calculate the 
induced current density, electric polarization, and magnetization as the expectation values of 
the corresponding operators. ) 



2.2 Microscopic response 

In order to determine the microscopic response, we only need to have the microscopic constitutive 
equation between /t and A, which is to be solved selfconsistently with eq. (fT3|) . 

The microscopic calculation of induced current density can be found in various textbooks as 
e.g., [ini [m [13 US], but the outline of the lowest order perturbation calculation is reproduced 
in Appendix B for the convenience of readers. 

The uj Fourier component of the induced current density is given in eq. (|65p . or in an integral 
form as 

/(r, Lo) = J dr'xir, r'; u) ■ A{r' , lo) (18) 
where the integral kernel is the nonlocal susceptibility of the form 

X(r, r'; cu) = - J] [5.(t^)i"o.(r-)/.o(r') + /i,(cu)/,o(r)/o.(r') ] (19) 

and gyiyi) and hyiyi) are defined in eq. (j66|) . The operator / is the A independent part of 
generalized current density ([59|) . 

The induced current density (165p is a linear combination of the factors {F^jy}, defined in 
eq. ([63]) . and the factors {F^^y) are determined by vector potential. The selfconsistent solution 
of the coupled equations for current density and vector potential goes as follows. 

Using eq. (l65p as the source term of eq. (ll3p . we can solve the M-eq. and obtain A(r,a;) as a 
linear combination of {FyQ^FQy\. Substituting this solution in the r.h.s. of eq. (j63p . we obtain a 
set of simultaneous linear equations of {F^q^ Fq^)^ the solution of which determines selfconsistent 
A and /t (in the region inside the matter). Once we have the selfconsistent A, we can use it 
to calculate /l via eq. (f65]) . and A at any point, e.g., response field outside a sample. In this 
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way the simultaneous equations of {F^,^} lead to all the information about the induced field and 
current density. This is the main framework of microscopic nonlocal response theory mentioned 
in Sec.l. 

The nonlocality of the susceptibility is the characteristic feature of this theory. Within the 
extension of the relevant wave functions, an EM field applied at a microscopic point r' can 
induce current density at different position r. Mathematically, the separable character of the 
integral kernel (susceptibility) plays an essential role to make this nonlocal scheme feasible in 
various manipulations, including the process of LWA to be discussed below. 

Another distinction from the macroscopic M-eqs. is the boundary conditions of EM field. 
To obtain the EM response in this microscopic scheme of nonlocal response, we do not need to 
consider the boundary conditions for EM field. (Boundary conditions are already used to solve 
the eigenstates of matter, which determine the nonlocal susceptibility.) Unique solution of the 
equations for {Fi^^} is obtained only by requesting a given initial condition of EM field, i.e., 
incident field. This way of solution is in sharp contrast with the conventional macroscopic M- 
eqs., which need the boundary conditions for EM fields. The various new features in microscopic 
response theory, suitable for the study of nanostructures in particular, are discussed in |18j . 

2.3 Long wavelength approximation of microscopic response 

As mentioned in the introduction, the logically straightforward, reliable way of macroscopic 
averaging is to apply LWA to the microscopic response described in the previous subsection. 
The validity condition of LWA depends on the system in consideration, so that it must be 
checked independently for each system. (See Sec. 3. 5 for more details.) In this section, we 
discuss the genaral form of the macroscopic M-eqs., assuming LWA as a good approximation. 

If LWA is valid, the spatial variation of vector potential A(r, uj) and induced current density 
I{r,io) will be weak in comparison with that of the matrix elements of the current density. 
Thus, the variables A and / are represented by their long wavelength components alone. The 
form of M-eqs. for A, eq. ()13p . is kept unaltered under LWA. In the Fourier representation, we 
have ^ 

{-^ + k^)Aik,u;) = —iTik,u;) (20) 

with the understanding that only small k components have appreciable amplitudes. A similar 
expression holds for the constitutive equation ([18]) as 

i{k,Lo) = -Y,\M^)iou{k)F^o{Lo) + K{uj)i^o{k)Fo^{u;)\ . (21) 

Here also, only small k components are considered to have appreciable amplitudes. 
The factor F^y{oj) can be rewritten as 

F^,{Lo) = Y,i^,{-k')-A{k\u:) , (22) 
k' 

which in general contains all the fc'-components, but, if LWA is valid, only small k's make the 
central contribution. 

For small fe, we may take the first few terms of Taylor expansion (around r = r) as 

V(fc) = Y^j dr exp[-^fc ■ r] I ^,{r) = ^""P^"^^ ' ""^ (V " " Q^.u) (23) 
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where = or z/ = 0, T4i is the normahzation volume, and 

V = /dr V(r) , Q^,u = J dr (r - r) V(r) (24) 

represent the moments of the electric dipole (El) and magnetic dipole (Ml) (plus electric 
quadrupole (E2)) transitions, respectively. We choose r for each transition " fi <-> z^" in such a 
way (e.g., at the center of impurity atom) that the moment Q represents its physical meaning 
correctly. Then, we obtain the relation between I{k, tu) and A{k', uj) as 

iik,u;) = ^ e^(fc'-^>^[5.(u;)(/o.-ifc-0o.)a.o + ifc'-0.o) 

+ K{uj){i^o - ik ■ Qm){Iou + ik' ■ Qou)] ■ ^(fe', w) . (25) 

This expression of the induced current density in LWA still allows the mixing of different wave 
vector components, which corresponds to a macroscopic body without translational symmetry. If 
the macroscopic medium obtained by LWA has a translational symmetry, as usually anticipated, 
we can keep only the k' = k term in the above summation. In this case, we have 

I(fc,w) = Xem(fc,w) • A(fc,u;) (26) 

where 

Xera{k,Uj) = — [gy{uj){I o„ - ik ■ Qq^){I ^ + ik ■ Q ^q) 

+ K(u)){IyQ - ik ■ Q^o){iou + ik ■ Qo^)] . (27) 

Here, we have replaced the factor l/V^ with the number density N^, of localized states (of 
impurities, defects, etc.) corresponding to the transition <-> z^, and the summation over v is to 
be taken only once for the same localised transitions at different sites. 

This is the general susceptibility of the present macroscopic scheme, and is the only suscep- 
tibility required to determine the complete (linear) response. Though this expression contains 
the both T and L components of /, only the T components are required for the selfconsistent 
solution with eq.(|20p. It contains the contributions from both electric (El) and magnetic dipole 
(Ml) transitions, together with their mixing terms. It should be noted that the mixing terms 
remain nonvanishing in the case of chiral symmetry, where each excited state > is active to 
both El and Ml transitions. In this case, we cannot properly define Xe and Xm, or e and fi. 
In the absence of chiral symmetry, on the other hand. El and Ml characters are not mixed, so 
that the k-linear terms vanish, and the susceptibility turns out to be a sum of El and Ml types 
of terms. This is the situation where we can use e and /x, and only in this case, the dispersion 
equation of the present formulation coincides with the conventional one in terms of e and fi, as 
will be discussed below in more detail. 



2.4 Dispersion equation of plane waves 

The macroscopic constitutive equation obtained above by LWA is described by the susceptibility 
tensor Xem{k,uj). Substituting the expression of I{k,uj) in the source term {Att/c)! of eq. ([20]) . 
we get the equation 

/ c^/c^ \ ~ Aire x-^, X T / N 
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where ^, rj are the two Cartesian coordinate axes perpendicular to k. The condition for the finite 
ampUtude solution is the vanishing of the determinant of the coefficient (2 x 2) matrix, i.e., 



det 



c2A;2 47rc 

— 5 1 5" Xcm[k,UJ) 



. 



(29) 



This is the dispersion relation in the present scheme of macroscopic M-eqs. It should be compared 
with the well known form of the dispersion relation in the traditional M-eqs. 



det 



. 



(30) 



Here also the T components of the tensors e and fj, should be inserted. 

Apparently, the two dispersion equations are different, because the contributions of electric 
and magnetic polarizations appear as a product in (j30p . while in (I29p as a sum (including an 
interference term). Moreover, the new result claims only one susceptibility, while there are two 
of them in the conventional formula. In view of the possible mixing of El and Ml transitions 
in the case of chiral symmetry, e and fi can have common poles, which leads to an unphyisical 
situation, i.e., the occurrence of second order poles in the product efi in spite of the linear 
respose. 

In the absence of chiral symmetry, however. El and Ml (+ E2) transitions are grouped into 
different excited states. Namely, there is no excited state > making both of I^o and Q^^q 
nonzero. In this case, we may divide the u summation into two groups, so that we have 



Xe 



Xe + Xr 



(31) 



where Xe and Xm are the partial summations over u for El and Ml (+ E2) transitions, respec- 
tively, defined as 



1 



Xe 



— 5] iV, [5,(a;)(fc • Oo.) • Quo) + ^v(^)(fc • O.o) " OoJ] 



(32) 
(33) 



for a unit vector k = k/\k\. In this case, the dispersion equation takes the form 



det 



(1 + AttXc + 4^x r 



, 



which should be compared with the traditional form of dispersion equation, ([30 

c'k^ 



det 



(l+47rxe)(l+47rx„ 







(34) 



(35) 



It appears that, even in this simplified case of the present formulation, the traditional form of 
dispersion equation cannot be recovered. However, if we use, instead of Xm, the more funda- 
mental magnetic susceptibility xb defined by M = xbB (see the next section, Sec. 3.1, for more 
details), we have = 1 + 47rxm = (1 — ^'^Xb)~^- Then, eq. ([35|) can be rewritten as 



det 



c^k^ 



(1 + 47rXe + ^TT—^XB, 







(36) 
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where the contributions of El and Ml transitions appear as a sum and the Ml term contains 
the factor of 0{k'^), as in eq. ([M]) . The relationship between (xe^Xm) and {Xe,Xm) is obtained 
from the comparison of 0{k^) and 0{k'^) terms of the two dispersion equations as 

Xe = Xe , XB = {u/ckfXui- (37) 

In this way the equivalence of (j34p and (j35p is demonstrated in the absence of chiral symmetry. 

The above argument gives the answer to all the problems 1) - 4) in the introduction: 1) The 
electric and magnetic susceptibilities (xe and xb) correspond to the first and second order terms, 
respectively, of the LWA expansion of the microscopic susceptibility x{f^ '''') in the absence of 
chirality. Namely, they are two tensors derived from a single nonlocal susceptibility. 2), 3) Since 
we do not separate / into different components, there arises no problem of non-uniqueness. 4) 
The (apparently) different wave number dependence of the two Ml transitions arises from the 
different means of description, i.e., Xm (or xb) in the case of spin resonance [8] and Xm in the 
case of orbital Ml transition [9]. 



3 Discussions 

3.1 Microscopic derivation of Xe and Xm 

To calculate induced electric and magnetic polarizations separately, we need the operator forms 
of these physical quantities. It is known that the decomposition of J into the current densities 
induced by electric and magnetic polarizations in the form of eq. (j52|) is possible also in operator 
form. If we use 

P(r) = / du ^eiri 5{r - uri) , (38) 
Jo ^ 

M{r) = u du e^ri x V£ 5{r — uri) , (39) 

Jo Y 

we can prove V • P = —p and eq.([52|) for charge neutral systems (See Sec. IV. C of [2]). 

Using these operator forms, we can calculate the induced electric and magnetic polarizations 
microscopically along the line of Appendix B. Then, the application of LWA to the microscopic 
constitutive equations gives the macroscopic (local) susceptibilities. Such induced polaizations 
are given as functions of A, since the interaction Hi^t is given in terms of A, so that the 
susceptibilities do not correspond to Xc or Xm- This way of calculation leads to the equivalent 
result to that of Sec. 2, because the calculated P{r) and M{r) gives the induced current density 
according to eq. (f52]) . 

In order to introduce the variables of electric and magnetic field explicitly in the interaction 
Hamiltonian, the following transformation of Lagrangian is known to be useful. Making use of 
the fact that the addition of a total time derivative of arbitrary function (of time and position) 
does not affect the least action principle of Lagrangian, we add the following term [21 [19] 

1-1 p 

F(t) = dr P-A (40) 

dt c J 

to the Lagrangian ([1]). The combination of this term with Hi^t ([6|) leads to a new form of the 
interaction Hamiltonian as 

HU = -J dr{P-ET + M-B} , (41) 
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where we have used partial integration and = —(l/c){dA/dt), B = V x A . 

The interaction Hamiltonian H^^^ and the operator form of P and M appear to be appropri- 
ate for the calculation of electric and magnetic susceptibilities. However, electric and magnetic 
fields cannot be applied independently for finite frequencies, and in chiral systems both of P 
and M can be induced by both E and B. This is not a situation to define electric and magnetic 
susceptibilities. 

Thus, the new form of interaction Hamiltonian does not guarantee the possibility of defining 
electric and magnetic susceptibilities in general. In the absence of chirality, however, there is no 
mixing of El and Ml transitions, so that P is induced by E alone, and M by B alone, which 
allows the definition of electric and magnetic susceptibility. In this case, however, a perturbation 
calculation similar to the one used for I(k,uj), (j26p . should give an induced magnetization in 
the form, not M = XmH, but M = xbE, since the interaction term (I4ip is linear in B. The 
poles of xb correspond, similarly as in the case of Xem, to the excitation energies of matter, i.e., 
magnetic excitations in this case. Then, the usual definition B = H + AttM and M = XmH 
lead to Xm = Xb(1 — 47rxB)~^ and /U = (1 — 47rxB)~"'^- The last relation has provided an essential 
key to prove the equivalence of eq. (p4l) and eq. (f35]) . The use of xb is a logical consequence, but 
there are so many literatures using Xm on behalf of xb- The reason for it is difficult to find out 
from the literature of M-eqs. However, this paper is not the first one to claim the use of xb- 
See [IT] for example. 

3.2 Comparison with other schemes of macroscopic M-eqs. 

There have been several proposals of single susceptibility schemes of macroscopic M-eqs. [10 ^ 111 1 
112] , which are apparently motivated by the argument of Landau-Lifshitz [4j, that, as frequency 
increases, magnetization tends to lose its physical meaning, i.e., it becomes meaningless to 
separate current density into the contributions of P and M, as in eq. (l52p . This is very similar 
to our motivation, but our statement is one step stronger than theirs. Namely, "for any frequency 
region, we should need only one vector field / as in microscopic M-eqs., and it should be possible 
to write macroscopic M-eqs. without separating / into the contributions of P and M" . 

In both of the proposals, they derive the microscopic forms of induced current density, from 
which one could derives conductivity and translate it into dielectric function. The calculation 
of conductivity is essentially equivalent to the content of Appendix B of this paper. As to the 
translation of conductivity into dielectric function, it is a standard matter for the case of achiral 
symmetry, but it is not clear if one can do it also for chiral case. The authors of [llj call their 
scheme ^^E,D,B" approach. In the same manner, we may call ours "(macroscopic) E,B" ap- 
proach. In [i2j, they describe the case of dipole approximation, which corresponds to our 0{k^) 
term of Xem- Our results in Sec. 2 show the importance of the whole set of "0{k^), 0{k^), 0{k'^)" 
terms of the Taylor expansion to obtain the consistent picture of single suscptibility scheme of 
macroscopic M-eqs., which is free from the incompleteness problems 1) - 4) in Sec.l. Both of 
[lOj and |12j stress the distinction between the presence and absence of spatial dispersion within 
the macroscopic M-eqs., but we understand this simply as the different order terms of Taylor 
expansion within LWA. In spite of the similarity in motivation and a part of formulas, neither 
of them show any systematic expansion of the microscopic constitutive equation to obtain a 
generalized form of macroscopic M-eqs. In this sense, we could say that the present method and 
result reaches a more complete level of single susceptibility scheme. 

There are a few points which hinder the proper evaluation of these proposals. In |llj . 
they claim the existence of "one-to-one correspondence" between their scheme and the usual 
{E, B, D, H} framework. If this means the equivalence of the two schemes, their result is 
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definitely different from ours. Our result is, not equivalent to, but more general than {E, B, 
D, H} framework, based on the recognition of the latter's incompleteness problems and of the 
necessity to rationalize them. In [12] , the microscopic matter Hamiltonian is described by the 
vector potential A defined as V x A = H. This should lead to the magnetic susceptibility Xm 
rather than to XB, yet "their magnetic susceptibility contributes to the single susceptibility a 
term of the type (c^fc^/a;^)xm" (p-92, [12]) as in eg. ([36]). The formulas related with magnetic 
current density in Sec. 1.5 of }12j 



is also difficult to understand. It is not clear if their claim is equivalent to ours as to the necessity 
to prefer xb over Xm- 

There is a different approach to macroscopic M-eqs. by Nelson [20], where he applies LWA 
to the Lagrangian of matter-EM field system, rewriting it into a "continuum" Lagrangian. 
The explicit use of LWA in the mathematical treatment is a common feature to the present 
theory, but the physical meaning is quite different. By the application of LWA to Lagrangian, 
the dynamics of matter is described only by the long wavelength (LW) components, i.e., the 
LW eigen modes of matter such as acoustic and optical phonons and excitons. Thus the only 
contribution to susceptibility is made from the LW modes of matter, i.e. the susceptibility 
has poles only at the frequencies of these LW modes. Since all the dynamical variables of 
short wavelength components are eliminated by the LWA of Lagrangian, there is no chance 
for localized eigen states of the matter to contribute to susceptibility. When we consider a 
problem, for example, of changing the refractive index of a material by adding impurities, the 
main change is caused by the localized excitations at the impurities. But they cannot be taken 
into account in Nelson's treatment, because they are not LW modes of matter. In contrast, our 
approach based on the LWA of microscopic constitutive equation takes all the contributions of 
the eigen modes of matter according to their weights in LWA, i.e., oscillator strengths. Thus, our 
approach provides continuous relationship between microscopic and macroscopic descriptions of 
EM response, including the method to evaluate the validity condition of LWA (see Sec. 3. 5). 

3.3 Comparison with Drude-Born-Fedorov equation 

Chirality or gyrotropy of a system appears as the rotation of polarization plane [6] and as the 
mixing effect between electric and magnetic dipole characters, e.g., of the excitons in CdS pTj 
and the coupled "Landau level - spin flip" transitions in GaAs [22J, and the Jones effect in 
atomic spectroscopy [23]. It arises from the interference of electric and magnetic polarizations. 
The description could be either microscopic [21], [22], [23], or macroscopic depending on the 
problems. The microscopic one is equivalent to the microscopic nonlocal response in Sec. 2. 2, 
and as macroscopic one we have two schemes, the present one in Sec. 2. 3 and 2.4, and that of 
DBF eqs. [6j. 

In this subsection, we want to show that the two macroscopic schemes cannot be equivalent. 
For this purpose, we use the dispersion relations. In our scheme, the dispersion relation is given 
as eq. (j29p . where the terms in the determinant other than c^fc^/w^ is a sum of single poles except 
for u = 0. This is a general aspect of linear response. The dispersion equation for DBF eqs. is 
obtained by solving the M-eqs. cV x H = dD/dt, cV x E = —dB/dt and DBF constitutive 
eqs. (j54p . The condition for the existence of finite amplitude solution of the coupled linear 
equations for E and H gives the dispersion relation. In the simplified case where e, /i, (3 are 
scalars (with uj dependece), the dispersion equation is given as 



J m 



cVxM = cVx Xm{H - A-kM) ~ C V X Xm-H" • 



(42) 




-2 



(43) 
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This is obviously different from the equation ()29p with respect to their pole structure on the r.h.s. 
Since El and Ml characters of matter excitations are mixed in chiral materials, the successful 
trick in Sec. 2. 4 to divide Xem into Xe and xb for achiral materials fails in this case. It is hopeless 
to rewrite Xem as a sum of El, Ml, and chiral components in accordance with DBF eqs., and to 
make the r.h.s. of eq. (l43p a sum of single poles. 

This means that the DBF eqs. cannot be justified from a microscopic basis. Though the 
parameter (3 can qualitatively describes the different phase velocities of right and left circularly 
polarized lights, we cannot obtain its quantum mechanical expression consistent with Xem- Thus, 
the use of the k-linear term of Xcm is preferable to take care of the chirality, rather than the 
DBF eqs. 

3.4 Influence on the studies of Left-Handed Materials (LHM) 

According to the results of the previous section, Veselago's definition of LHM (e < 0, ;U < 0) 
|16j needs to be revised because of the limited condition to allow the use of e and /i. A more 
general definition would be "the occurrence of a dispersion branch with x < 0", where 
t;ph and Vg are phase and group velocities, respectively. The part of positive Vg in the negative 
k region is important because the plane wave on this part of the branch should be connected to 
the incident EM field via the boundary conditions at the interface. When two (or more) plane 
waves exist in the medium at the frequency of an incident wave, the general criterion for the 
correct choice of the branch is that the connected wave in the medium should not diverge as 
the size of matter becomes larger. From a general consideration of the dispersion curves in the 
complex {uj,k) plane \2^, it is concluded that the choice of the branch with positive Vg meets 
this criterion. Physically, this means a simple fact that the plane wave with postive Vg deays in 
the direction of its propagation. Therefore, it applies to any dispersion curves of either right- or 
left-handed character. 

LHM can be made of either chiral or achiral materials. Though it is recommended to use Xem 
in general, it is possible to use the conventional "e,//" scheme for achiral materials. In this case, 
however, it is also required to assign the magnetic transition energies to the poles of, not Xm, 
but Xb- Though the assignment to xb is logically correct, that to Xm has been very frequently 
made in textbooks and papers. In order to show the possibility of a serious difference due to 
this distinction, we give a measurable example below. 

A typical LHM behavior is expected in an achiral material when a magnetic transition occurs 
in the frequency region where 1 -|- (47rc/u;^)xc (= ^b) < . The two different assignments are 

(A) Xm = (44) 

(B) XB = (45) 

where b and are the strength and resonant frequency of the magnetic transition, respectively. 
The corresponding dispersion relations are 

A = eb 1 + — , 46 



B = eb (1 — )-^ . (47) 



The solution of this equation gives the dispersion relation k = k{uj) and refractive index n = 
ck/uj. For the positive Vg branch of LHM behavior, k and n are negative. The reflection 
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Figure 1: Dispersion curves with fgfph < 0. The parameter values eb = — 1 and Anb/uJo = 0.001 
are used. The abscissa (k) and ordinate (a;) are normahzed hy uj/c and o^o, respectively. Two 
curves correspond to the cases (A) and (B) in the text. The forms of the curves are quite similar, 
except for their positions with respect to the resonant frequency luq. 
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Figure 2: Reflectance spectra of the two cases in Fig.l, showing transmission windows in the 
total reflection range. Two curves are almost identical, except for their positions with respect 
to the resonant frequency lo/loq = 1. 

coefficient of a semi-infinite system for normal incidence is given as R = \{n + l)/(n — 1)^ for a 
medium of LHM, with n < 0. 

Dispersion curve and reflectivity spectrum are shown in Fig.l and Fig.2, respectively. We 
have a branch with VgVph < as expected, and this branch opens a window (a dip) in the 
total reflection spectrum (R = 1). Since the value of n along the dispersion curve changes 
continuously between and — oo, the dip of reflectivity can become as deep as (n = — 1). The 
remarkable point is the relative position of loq with respect to the reflectivity dip. In (A), loq 
occurs on the lower frequency end of the window, and in (B) on the higher frequency end. This 
could be checked experimentally, via an appropriate model system, e.g., a magnetic resonance 
of well-defined impurities overlapping with an El type phonon resonance. 

3.5 Validity of LWA 

In our derivation of the macroscopic M-eqs., we have just assumed the validity of LWA for the 
system in consideration, and applied LWA to the microscopic response. The validity condition 
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of LWA is not provided by the macroscopic theory itself, but must be checked independently 
through the microscopic scheme described in Sec. 2. 2, where all the wavelength components are 
included in the response field. Only when the amplitudes of the short wavelength components 
are small enough in comperison with the LW component, LWA is a good approximation. 

The validity condition of LWA is closely connected with the resonant or non-resonant con- 
dition of the optical process in question. In resonant optical processes, only a few excited states 
will be resonant to incident frequency, and they will make the main contribution to the induced 
current density. Thus, the current density acquires a characteristic microscopic spatial structure 
reflecting the quantum mechanical wave functions of the resonant excited states. All the other 
non-resonant states make the contributions of more or less comparable amplitudes, so that their 
superposition will have no particular microscopic structure. Therefore, non-resonant processes 
could generally be handled by macroscopic scheme, and resonant processes should mostly be 
treated by microscopic theory. However, the resonant processes due to uniformly distributed 
well-localized impurity or defect states may well be treated by the macroscopic theory, as long 
as one considers the average properties of these localized states. 

In the case of resonant optical processes in confined systems, the breakdown of LWA can 
occur rather quickly as we increase the size of confinement starting from an atomic scale. (See 
sec. 4.1 of |18j.) In this case, one should compare the confinement size, not with the wavelength 
of the resonant light in vacuum, but with that in matter which is shortened by the factor of 
background refractive index. 

3.6 Simplification of Xcm with parameters 

The merit of the conventional macroscopic M-eqs. lies in the description of matter in terms of a 
few material parameters. Namely, for each material, one assumes certain simplified forms of e, 
/i, and chiral admittance containing some parameters. Thereby one could consider resonant as 
well as nonresonant behavior of matter, so that one can handle a wide variety of EM responses. 
These parameter values can be determined from the comparison between the solution of the 
M-eqs. and certain appropriate experiments. 

In this subsection, we discuss how far this kind of approach can be justifed from the results 
of this work. While the conventional approach uses three susceptibilities e, /x, and chiral ad- 
mittance (3 defined for different pairs of vector fields, our new result requires the use of a single 
susceptibility Xem between I and A, which contains the three tensor terms of the order 0{k^), 
0{k^) and 0(A;2) as 

Xem = Xel + ik Xchir + Xml ■ (48) 

Each term on the r.h.s., to be obtained from eq. (l27p . consists of a summation over excited 
states u, which could be simplified as a constant (plus a few resonant terms). This will reduce 
the number of freedom, facilitating the analysis of experimental results. Whether or not the 
parameters in the simplified expressions can be treated as "free" should be judged from the 
original expressions. In chiral materials, for example. El and Ml transitions mix each other, 
so that all the three terms have common single poles. Therefore, if one keeps a pole in the 
simplification process, one should keep it for all of Xei, Xchir and Xmi- 

It would generally be safe to use Xei, Xchir and Xmi with common resonances. Especially, 
the models of LHM without inversion symmetry such as split ring resonators would generally 
need such a treatment. Only in nonchiral materials we can use "e and /i" with free parameters. 
But we should still be careful in ascribing a magnetic transition energy, not to the pole, but to 
the zero of /i. 
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4 Summary 



The main points of our logical steps and important results are as follows. 
Logical steps 

1) The derivation of the conventional macroscopic M-eqs. from microscopic basis is not 
complete because of the lack of uniqueness in dividing current density into the components 
arising from electric and magnetic polarizations, and of the apparent inconsistency with 
microscopic response about the number of constitutive equations. 

2) As a new method of macroscopic averaging, we apply LWA to the microscopic nonlocal 
constitutive equation (and the microscopic M-eqs.) without any other assumption. 

3) The whole derivation is made for the precisely defined Hamiltonian of charged particles 
interacting with oscillating EM field, where an explicit consideration is added about spins 
in the Hamiltonian and the current density operator. 

Important results 

a) The macroscopic M-eqs. obtained from the above mentioned procedure retain the same 
form as the microscopic M-eqs., written in terms of E and B (or A and 0), which require 
only one susceptibility tensor Xem- 

b) This susceptibility tensor contains the contribution from both elecrtic and magnetic po- 
larizations, together with their mutual interference, which is in generl not separable into 
different components of conventional type. 

c) Since it is not necessary to introduce P and M, this scheme is free from the problems 
about uniqueness and consistency inherent to the conventional one. 

d) Matter with chiral symmetry, where El and Ml transitions are mixed, should be treated by 
the present framework. The phenomenological scheme with additional chiral admittance 
of Drudc-Born-Fcdorov cq. cannot be reproduced from the microscopic theory. 

e) The dispersion equation is det|(cfc/a;)^ — 1 — (47rc/ti;^)xem| = in general. This can be 
reduced to the well known form det| (c/c/w)^ — €//| = in the case of achiral symmetry. 
Thereby, it is essential to use xb (defined by M" = xb-B) and /x = (1 — 47rxB)~^- 

f) The linear response coefHcient of the magnetic transitions in nonchiral matter is not Xm 
(defined by M' = XmH) but xb- An observable difference in these assignments is shown in 
connection with LHM behavior of a magnetic transition in the frequency region of e < 0. 

g) The susceptibilities Xc and xb correspond to the first and second order terms, respectively, 
of the LWA expansion of the microscopic susceptibility x(^)^')'^) in achiral materials. 

h) Due to the restricted condition to use e and /j,, a better definition of LHM would be the 
occurrence of dispersion branch with " (phase velocity) x (group velocity) < 0" . 

i) For the study of LHM, it is necessary to check whether or not a proposed model structure 
allows the use of e and /i, and to ascribe the resonant frequncy of magnetic excitation, not 
to Xmj but to Xb- More generally, it is recommended to use Xem- 
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Appendix A. Fundamentals of M-eqs. 

Though the contents of this appendix is well-known, we give them explicitly in order to 
provide a basis to check all the formulas of this paper and to fix notations and units. 

Microscopic M-eqs. are the equations for two field variables, electric field E and magnetic 
field (flux density) B as 

V-E = A-Kp, V-B = 0, V X B = —J + -^, VxE = --^. (49) 

c c at c at 

The source terms p and J are microscopic charge and current densities, respectively, satisfying 
the continuity equation, i.e., charge conservation law, 

V-J + |^ = 0. (50) 

These equations determine E and B for given p and J. To describe the EM response of 
matter, we need additional relationship to express p and J in terms of E and B. This kind of 
relationsip is called constitutive equation, and, in the present case, we generally need only one 
integral equation relating J and or S". Since p and J are related via continuity equation, 
and E and B via Faraday's law, only one relationship between two vector fields is sufficient. 
Thus, the number of required susceptibility tensor (integral kernel) is one. 

Macroscopic M-eqs. are used for the description of EM field in macroscopic matter, i.e., gases, 
liquids and solids, so that they employ macroscopic variables (usually electric and magnetic 
polarizations plus macroscopic current density) rather than the microscopic charge and current 
densities mentioned above as source terms of the equations. The central idea to derive the 
macroscopic from microscopic M-eqs. is the "macroscopic averaging of microscopic variables", 
which is made "over a volume larger than atomic scale but smaller than light wavelength" to 
extracts the macroscopic components of field and matter variables. 

The standard way to derive macroscopic M-eqs. from microscopic M-eqs. is to separate 
certain parts of charge and current densities as the contributions from electric and magneteic 
polarizations as 

P = Pt + Pp , Pp = -V • P , (51) 

dP 

J = j^ + cV X M + — , (52) 

where pp represents polarization charge density, c V x M and dP/dt the current densities due 
to magnetic and electric polarizations, respectively. The polarization charge density is defined 
for the part of neutral charge distribution, i.e, J dr pp = 0. The remaining (non-neutral) part 
of the charge density contributes to pf The true charge density pi and the current density Jc 
caused by its motion satisfy the continuity equation of the form of eq. (|50p . In terms of the 
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electric polarization P and magnetic polarization M, we introduce the new field variables D 
and H as D = E + At:P and H = B — AttM. Then, the microscopic M-eqs. are rewritten as 



V-D = A7Tpt, V-B = 0, Vxi/ = — Jc + - — , VxE = —— (53) 

c c at c ot 

This set of equations are the macroscopic M-eqs. where all the field variables and source terms 
are macroscopically averaged quantities consisting of long wavelength components alone. 

The constitutive equations for macrsocopic M-eqs. are required for both P and M as 
P = XeE, M = XraH in the case of linear response. The material parameters Xe and Xm are 
called electric and magnetic susceptibility tensors, respectively, and they define the dielectric 
constant and magnetic permeability as e = 1 + 47rxe, /i = 1 + 47rxm [D = eE, B = fiH). 

There are low symmetry systems where electric and magnetic polarizations are mixed with 
one another, and the plane of polarization of a linearly polarized light rotates as it propagates. 
This can occur either naturally or by the application of static electric and/or magnetic field. 
This kind of system is called chiral, optically active and/or gyrotropic, and the phenomenological 
way of its description is known as Drude-Born-Fedorov (DBF) constitutive equations [6], [7], 
which takes into account an additional susceptibility called chiral admittance /3 as 

D = e{E + l3V X E), B = fi{H + /3V x H) . (54) 

In this case, we have two constitutive equations with three susceptibilities. 



Appendix B. Derivation of microscopic susceptibility 

Here we calculate the expectation value of current density operator with respect to the 
matter wave function at time t, < ^{t)\I{r)\^{t) >. The time evolution of ^ is governed by the 
Schrodinger eq. ifid^/dt = Hq + H^at, where the matter Hamiltonian Hq and the matter-EM 
field interaction H\at are defined in eqs. (llOp and (jllh . respectively. We neglect the 0{A^) term 
in the interaction Hamiltonian, since we are interested in linear response. 

Using the interaction representation ^ = exp{—iHoT/h)^, ihd^/dt = H'{t)^, where 

H'{t) = exp{iHoT/h) i/int exp{-iHoT/h) , (55) 
we can easily obtain the solution to the first order of A as 

§(t) = §(-cx)) - ^ /* dr H'{t)'^{-oo)+--- . (56) 

n J-oo 

The interaction i^int is adiabatically switched on in the remote past, where the matter state is 
assumed to be in the ground state |0 > of the Hamiltonian Hq^ i.e., ^'(— oo) = |0 >. 

The A linear terms of the expectation value < ^ {t)\I {r)\^ [t) > originate from those of ^{t) 
and of the operator /. The latter is in the velocity term in the orbital current density, eq.®. 
Since the velocity of each particle in the presence of vector potential is given as 

1 6/? 

vi = —[p^--A{r^,t)], (57) 
m£ c 

we can rewrite J as 
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where N{r) = ^^^(ef/m^) 6{r — Vf) . 

Let us separate the A dependent part of /, i.e., the second term on the r.h.s. of eg. ([58]) . as 

I = i --N{r)A{r,t) . (59) 
c 

Then, the A Hnear part of the induced current density < ^ {t)\I {r)\^ [t) > is given as 

I{r,t) = -- < 0\N{r)\0 > A{r,t) - ^ f dr < 0| [i{r)H'{T) - H'{T)I{r)] |0 > (60) 

Using the eigenstates {{i^ >, |/u > etc.} of Hq, we can evaluate the matrix elements of H'{t) as 

< fi\H'{T)\v >= -^exp[i{Ef, - E^)T/h] J dr < fi\I{r)\iy > -Air, t) . (61) 

This form, together with the adiabatic switching factor exp(0~''r) in -ffinti where 0"*" is an in- 
finitesimally small positive number, allows us to carry out the time integral in eq. ()60p . The a; 
Fourier component of i"(r, t) gives the microscopic constitutive equation 

/(r,u;) = --<0|iV(r)|0> A{r,uj) + -Y,[9v{^)Iov{r)FM{uj) + K{iv)I^Q{r)FQ^{uj)] , (62) 
where 

F^^{uo) = Jdr I^^{r) ■ A(r, w) (^ = 0, or = 0) , (63) 

and 

The first term on the r.h.s. of eq. (l62|) is the contribution of the ground state charge density. 
This term is important in several situations, such as X-ray diffraction, low frequency conductiv- 
ity, etc. For resonant optical processes in nanostructures which is the most important playground 
for the microscopic resonse theory, however, the main contribution comes from a few resonant 
excited states and the first term makes one of very many nonresonant contributions which is 
altogether treated as a background susceptibility. For the derivation of macroscopic M-eqs., 
which is the main purpose of this theory, it would be enough to retain the main part of the 
first term, (— e^no/moc) A, where e, niQ, and uq are the charge, free mass, and averaged ground 
state number density, respectively, of electrons in matter. In this case, it is possible to rewrite 
the first term into a similar form as the second one, which results in a simple expression of the 
current density 

I{r,Lo) = -J2[M^)iou{r)F^o{uj) + K{uj)i^o{r)Fo^{Lo)] , (65) 
where ^ ^ 

giy{^) = 9y{^) - TT- , = - — , (66) 

and EyQ is the matter excitation energy from the ground to the i^-th excited state. The terms 
with the factor I/E^q are those arising from the first term eq. (|62p . For this manipulation, we 
make use of the commutation relation between current density and dipole density operators and 
LWA is used for vector potential. The details is given in Sec. 2. 4 and 2.5 of [18]. 

When matter state is described by ensemble more generally, e.g., a canonical ensemble at 
temperature T, one can use the perturbation expansion of the equation of motion of density 
matrix [12], [TTj. (The result eq. ()62p corresponds to the canonical ensemble of T = K.) Since 
the (r, r') dependence of the microscopic susceptibility is also separable in this case, the process 
of applying LWA is very similar to that described in Sec. 2. 3. 
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